Evolution of the Pseudogap in a polarized Fermi gas 
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We calculate the single particle spectral density of a normal (non-superfluid) two component gas 
of fermions in the BCS-BEC crossover within a T-matrix approximation. We review how non- 
condensed pairs lead to a spectral density reminiscent of the ordered state, and explore how a 
gap-like feature in the spectrum evolves as one changes the polarization of the gas. As the gas is 
polarized we find that this pseudogap becomes more diffuse and moves away from the Fermi level, 
reflecting the fact that fewer pairs are present but that they still play an important role in the 
excitations. 
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Two closely related themes in modern condensed mat- 
ter physics are how correlations manifest themselves in 
unordered states of matter, and the existence of materials 
with unordered ground states. These themes are ubiqui- 
tous, playing a role in studies of spin systems, transition 
metal oxides, and cold gases: they are key to phenom- 
ena such as spin liquids and resonating valence bonds 
[l[ . One important idea to emerge in this context is that 
of a "pseudogap": the fact that under appropriate cir- 
cumstances the normal state displays a suppression of 
the single particle spectral density near the Fermi level, 
reminiscent of the gaps seen in ordered states of mat- 
ter. Here we present a theoretical study of the pseudo- 
gap in a strongly interacting Fermi gas, such as 6 Li or 
40 K atoms. In particular, we show how this pseudogap 
structure evolves as one polarizes the gas. 

Over the past decade, the pseudogap system which has 
received the most attention has been high temperature 
cuprate superconductors Q . In that system the origin of 
the pseudogap is widely debated, with some arguing that 
it is actually a real gap associated with non-superfluid or- 
dering. Less controversially, there are at least two other 
systems where vestiges of ordering dominate the prop- 
erties of the normal state: atoms in the "BCS-BEC" 
crossover [3rfl5l [l7l - r20j , and electrons in one-dimensional 
charge-density wave materials [2l|. The former system 
is ideal for studying pseudogap phenomena: one has a 
wealth of "knobs" for adjusting the parameters, and a 
number of very powerful probes. For example, we will 
calculate the single particle spectral density, which can 
be extracted from RF spectroscopy experiments fl8l - l2fjj . 

These experiments on cold fermi gases have dra- 
matically shifted our perspective on superconductiv- 
ity/superfluidity [20, 22]. They have confirmed the ideas, 
largely developed in the 1980's, that the superfluidity 
seen in 4 He - where bosonic atoms condensed - is contin- 
uously connected to the superfluidity seen in 3 He - where 
superfluidity results from a fermi surface instability [23l — 
[251 ] . In the cold gas experiments a magnetic field is used 
to tune the energy of a two-body state [26j : taking it from 
a negative value (where it is a bound state) to a positive 
value (where it is a scattering resonance). In the former 
"BEC" limit, the low temperature system is described as 



a gas of bosonic dimers, while in the latter "BCS" limit 
it is a gas of fermions, interacting via an effective attrac- 
tive interaction. The superfluid ground state smoothly 
evolves as one changes the magnetic field: the molecu- 
lar pairs on the BEC side simply grow into large Cooper 
pairs as one approaches the BCS limit. In this crossover, 
one of the most interesting point is where the free-space 
bound state has exactly zero energy. Here interactions 
provide no length-scale, and the scattering cross-section 
saturates the bounds set by the unitarity of the S-matrix. 
At this "unitary point," the thermodynamics potentials 
are "universal" functions of the dimensionless combina- 
tions A 3 n,j, where n a is the density of atoms with spin a 
and A 2 = 2irh 2 /mkFjT is the thermal wavelength, where 
h = h/2n is the reduced Planck's constant, m is the 
atomic mass, fee is Boltzmann's constant, and T is the 
temperature [27j . The experimental systems are remark- 
ably stable near unitarity. 

One difficulty in studying the pseudogap is the ab- 
sence of a single universally accepted definition. In the 
context of high temperature superconductors the term is 
typically used to refer to a number of phenomena in the 
underdoped normal state. Only some of these proper- 
ties are found in the BCS-BEC crossover (see [lg, [l7[ for 
reviews) . Throughout this paper I will define the pseudo- 
gap in terms of a dip in the density of states. This tends 
to be one of the weaker definition of the pseudogap in 
the BCS-BEC crossover, and the one which is common to 
all sufficiently sophisticated theories of the normal state 
near T c . Other possible definitions, such as bimodality 
of the momentum resolved single particle spectral den- 
sity at the Fermi energy @-|a5, or "back-bending" of the 
peak in the single particle spectral density [HI , are more 
stringent: not all theories display them. For example, 
the T-matrix theory used in this paper does not show a 
definitive "back-bending" structure Q, especially when 
the gas is polarized. In response to this plurality of def- 
initions, Tsuchiya et al. [5[ introduced the phenomono- 
logical concept of two "pseudogap temperatures" , one as- 
sociated with the appearance of a dip in density of states, 
and another with the appearance of "back-bending" . 

Alternatively, some use non-spectroscopic signatures 
to define the pseudogap. For example, Shin defines the 
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FIG. 1: Pseudogap in the density of states of the (two com- 
ponent spin balanced) normal state just above T c , calculated 
within our T-matrix approximation. Vertical Axis: density 
of states p(E) measured in terms of the density of states at 
the Fermi surface of a non-interacting gas po(Ef)- From top 
to bottom (at E = 0), l/k F a = -2.6, -0.7, -0.3,0.0. All en- 
ergies are measured from the chemical potential. The partial 
self-consistency enhances the suppression of spectral weight 
for strong interactions, broadening the pseudogap (see related 
figures in Refs. [ISiGl, which show how varying degrees 
of self-consistency influence the spectra). 



pseudogap state in terms of the equation of state [29j . 
According to his definition, the temperature dependence 
of the pressure found by Nascimbene et al [30j is incom- 
patible with a pseudogap. This interpretation is quite 
contentious: both Chien and Levin [28[ and Perali et al. 
[l8| point out that strong pairing correlations in the nor- 
mal state can lead to gap-like features in the spectrum 
yet produce thermodynamic functions which are consis- 
tent with the experiments. 

The experimental evidence for the pseudogap in Fermi 
gases comes from RF spectroscopy measurements fl8j . 
In those measurements one extracts the trap averaged, 
but momentum resolved, density of states. At each k, 
the experimentalists find the energy, Eh, for which the 
density of states is maximal. In both the superfluid and 
normal state they find this dispersion is non-monotonic. 
Fitting the dispersion to a BCS form E 2 = (k 2 — fi) 2 + 
Ap g , with free parameters /i and A pg , gives an estimate 
of the "pseudogap" energy scale, A pg . 

The physical picture used to understand the pseudogap 
in the BEC-BCS crossover is that of "preformed pairs" . 
This is conceptually clearest deep on the BEC side of 
resonance, where the normal state is described as a non- 
condensed gas of bosonic dimers (the "preformed pairs" ) . 
Any single-fermion excitations requires breaking a pair: 
either by supplying the energy with your probe, or by 
relying on a thermal fluctuation. The former excitations 
have gaps, while the latter ones are exponentially rare 
for T below the binding energy Et> = k^T* . Thus there 
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FIG. 2: Density of majority species (solid) and minority 
species (dashed) states in the zero temperature normal state 
of a spin imbalanced Fermi gas. Top: l/(kffa) = —0.7, Bot- 
tom: l/(kf[a) = —0.35; nj,/n-|- = 0.25. Energies are mea- 
sured from the chemical potential of each species, thin lines 
show noninteracting density of states. 



is an exponentially small density of states at the Fermi 
energy. As one approaches the BCS side of resonance, the 
temperature scale for pairing, T*, drops towards T c , and 
the pseudogap regime T c < T < T* becomes vanishingly 
small. 

Here we ask what happens to the pseudogap when the 
atomic gas is "spin imballanced" or polarized, meaning 
that > nj.. Experimentally one sees that as the polar- 
ization increases the transition temperature drops. This 
trend is not surprising: fewer 4,-particles means that it 
will be harder to form pairs. When n^/n^ = 7 exceeds 
a critical value 7* ~ 2.3, the gas remains normal even at 
zero temperature |30l - [32l |. 

Does this zero temperature normal gas have a pseu- 
dogap? Given that one typically interprets the pseudo- 
gap in terms of pairs, the presence of a pseudogap would 
nominally imply that there are non-condensed pairs in 
this zero temperature gas. A zero temperature non- 
condensed Bose gas is extremely exotic: Mott insulators, 
and bosonic fractional quantum Hall states, are two ex- 
amples that come to mind. 

Since such an exotic state seem unlikely, one expects 
not to see a pseudogap in the polarized gas. How does 
this happen? Where does the pseudogap go as one con- 
tinuously tunes the system from the high temperature 
unpolarized normal state to the low temperature polar- 
ized normal state? 

Our answers to these questions are summarized by the 
density of states in Figs. Q] and [5J Figure Q] illustrates 
that there is a significant depression in the unpolarized 
normal-state density of states at the Fermi energy when 
T > T c . This depression is how we define the pseudogap. 
Fig. [2] shows that in the zero-T spin- imbalanced gas this 
depression is pushed away from the Fermi energy. The 
spectral dip is at positive energy for the minority species 
density of states and negative energy for the majority 
species. These very broad depressions are continuously 
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FIG. 3: (Wide) Single particle spectral density in the normal state of an unpolarized two-component Fermi gas with T 
slightly above T c . Darker colors represent higher density of states. Left to right - weak to strong interactions 1/k/a = 
-2.6,-0.7,-0.3,0.0. 



connected to the spectral dips in Fig.[T] The remainder of 
this paper describes the physics behind this finite energy 
pseudogap, and details our calculation of the spectral 
density. We will argue that like the pseudogap in Fig. [TJ 
the spectral features in Fig. [5] arise from pairing physics. 
No pairs occur in the ground state, but the pairs play an 
important role in the excitations. 

We take a single-channel model, with Hamilto- 
nian H = K + V, where K = J^ka e kaa\ak and 
V = 9j2k P q a lt a ti a p-ql a k+Pl-> tne dispersion is t ka = 
k 2 /2m — [i a and the coupling constant g is related to 
the scattering length by g~ x — m/ (4nh 2 a s ) — yJ2 q f?- 
We use a T-matrix approximation to the single-particle 
self-energy, 



E t (fc,u;„) 
@fe(w n ) 



1,"n 



4irh 2 a 



(1) 



where is the volume of space, /3 = \/ksT is the inverse 
temperature, the frequency sums run over Matsubara fre- 
quencies, uj n = 2irim/kBT. The free Greens functions 
are 



(2) 



Physically this approximation amounts to solving the 
two-body problem, taking into account the presence of all 



of the other atoms only through Pauli blocking. This is 
the minimal model for capturing the BEC-BCS crossover: 
the corresponding free energy was introduced by Nozieres 
and Schmidt-Rink [25|. and it has been used/modified by 
numerous authors [44l3l [23l l32l l33j . Mathematically it 
can be derived in a functional integral formulation by tak- 
ing Gaussian pairing fluctuations about the noninteract- 
ing normal state. Although one has no a-priori reason to 
expect that this approximation is quantitatively accurate 
when the interactions are strong, it appears to work sur- 
prisingly well: it predicts T c « 0.2_Ey at unitarity, while 
Monte-Carlo [H finds T c = 0.152(7)^/. As described 
below, we will be using a modified version of this the- 
ory for which quantitative predictions are less accurate, 
but which avoids some unphysical features. For qualita- 
tive matters, such as the one concerning us, our modified 
theory should be reliable. The theory in ((T|) has been 
used extensively to investigate pseudogap physics in the 
unpolarized case [5| , where it displays features which are 
similar to the fully self-consistent T-matrix approxima- 
tion [Io|- As discussed by Chien, Guo, He and Levin 
different levels of self-consistency yeild stronger/weaker 
spectral features. The simple approach in Eq. ([I} yields 
a particularly difuse pseudogap, which is further broad- 
ened by our partial self-consistency. 

The use of free pr opegato rs in ([1]) leads to some un- 
physical results [32|, |3jj|, |36|. The key problem is that 
when one considers Pauli blocking in the two-particle 
collisions, one uses the densities of free particles with 
chemical potentials /i-j- and /ij,. These densities are much 
smaller than those of interacting particles at the same 
chemical potential. To mitigate this problem, we follow 
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FIG. 4: Single particle spectral density in the normal state of 
a T=0 two-component Fermi gas with l/kf^a — —0.035, and 
ni/rif — 0.25. Left: majority (t), Right: minority (J.). 

Schneider et al. [35[ and shift [i G in ([2]) so that = 

GS^E^Gf CO and n„ = (Z^)" 1 G fe K0 
are equal. We view this as a form of self-consistency 
[37j, where we approximate our self-energies as a con- 
stant. Perali et al. use a similar approach, but with a 
slightly different renormalization procedure fl^ . While 
it avoids some unphysical results, the quantitative pre- 
dictions of the partially self-consistent theory tend to be 
less accurate than the non-selfconsistent theory. For ex- 
ample, within our approximations we find T c = O.bEf. 
Except for some additional broadening, the rough quali- 
tative features of Figs. Q]|3] are the same with and without 
self-consistency. The spectral tails in Figs. [S][H] differ in 
the two theories. 

We tabulate the imaginary part of the self-energy, V = 
2Im(E), via 

/d 3 k 
-^—-^ A k+P (uj + efcf) [/ + g] , (3) 

where A k (uj) = 2Im(T fc (w)), / = l/(e /3efc t + i) ) an d g = 
]_^ e /3(w+cfct) — 1). We then use the Kramers-Kronig rela- 
tionship to extract the real part. From the single particle 
Greens function, G _1 = Gq 1 — E, we extract the spectral 
density Ak a {uS) = 2ImCr CT (fc, uj), which is the number of 
single-particle states which have energy uj and momen- 
tum k. The density of states is p a {E) — J2k Aka(E). 
Sample density of states arc shown in figure Q] for tem- 
peratures slightly above T c . As interactions are made 
stronger one sees a pronounced dip in the density of 
states near the Fermi surface. In figure |3] the full mo- 
mentum resolved spectral density is shown. As interac- 
tions grow, the bottom of the band is shifted downward 
- an effect analogous to the "Hartree" shift seen at weak 
coupling. The higher energy states have a smaller shift. 
The break-point, near the Fermi energy, has low spectral 
weight. If one looks carefully at these figures, especially 
the third from the left, one sees a faint downward dispers- 
ing band. These states can be interpreted as excitations 
where one creates a positive energy (ie. non-condensed) 
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FIG. 5: Partially self-consistent T- matrix calculation of T = 
occupation numbers for majority (solid) and minority 
(dashed) atoms as a function of momentum k scaled by the 
majority species fermi momentum kff. Dotted lines show 
the location of kf-f and kf±. Main figure shows occupation 
numbers scaled by (k/kf) 4 in order to illustrate their asymp- 
totic behavior, while the inset shows the unsealed occupation 
numbers. Here l/(fc/^a) = 0.35 and n^/rif = 0.25. 



"pair" . Crudely, one adds a f particle of momentum k to 
this band by adding a pair of momentum p (with energy 
-Epair(p)) and removing a J, particle of momentum p — k. 
This changes the system's energy by E va i T {p) — e p -ki- As- 
suming that the bulk of the spectral weight comes from 
the lowest energy pairs, which are at p — with vanish- 
ing energy as T — > T c (ie. impair (0) ~ 0), one should find 
single particle states near the line uj = — efcj.. In anal- 
ogy to what happens in the superfluid, the hybridization 
between these "hole" states and the regular "particle" 
states where u = e^f gives rise to the pseudogap. On 
the BEC side of resonance, which we do not discuss here, 
-Epair(O) < 0, and the entire downward dispersing band 
lies at negative energies. 

This "hole" branch is quite diffuse in the normal state. 
Its breadth comes from both the fact that one can add 
pairs with p ^ 0, and the fact that away from on the 
BEC side of resonance, the pairs do not have particularly 
sharp energies. This latter point is mathematically seen 
in the structure of the T-matrix. The T-matrix does not 
have a pole on the real axis, though near T c it becomes 
sharply peaked. In that regime our "cartoon" gives a 
reasonable picture of the underlying physics. We must 
add Schneider et al.'s caveat that one can get downward 
dispersing spectral weight merely from short range corre- 
lations [351 ] , though this argument mainly applies to large 
wave- vectors. On a purely technical note, our partial self- 
consistency makes these spectra more diffuse than those 
in @. 

As one polarizes the gas by increasing ^ and decreas- 
ing /ij,, the cartoon argument says that the crossing point 
moves: both the f-particle and 4,-hole branches move to 
lower energy, so the crossing moves to negative uj in the 
majority (f) spectral density . For the 4,-atoms the op- 
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FIG. 6: Partially self-consistent T- matrix calculation of T = 
RF spectrum / for majority (solid) and minority (dashed) 
atoms in the normal state of a polarized Fermi gas with 
l/(kffa) = 0.35 and n^/nf = 0.25. Spectrum is normal- 
ized by the peak value of the majority species intensity / max . 
Similar figures appear in [35| 



posite occurs, and the crossing in the minority spectral 
density is at positive u). This scenario resembles the 
"breached pair" states introduced by Liu and Wilczek 
[HI , but this is in the normal state. The physical picture 
is that there are initially no pairs, but one can create 
pairs while adding a particle. As discussed by Fumarola 
et al [39( , once you polarize the gas there is a sharp pole 
in the T-matrix (but the total weight in this peak is not 
necessarily large). 

Our calculation of the density of states is consistent 
with this "excited pairs" scenario. In Fig. [21 one sees that 
in the T = polarized normal state the majority species 
density of states has a slight dip at energies below the 
Fermi surface, and in the minority spectral density the 
dip occurs above the Fermi surface. These dips (particu- 
larly the one in the majority density of states) are much 
more subtle than those in Fig. [T] In the momentum re- 
solved spectra shown in Fig. |4l one sees that the dips are 
caused by the extremely short lifetime of the quasiparti- 
cle excitations in those energy ranges. This behavior is 
what one would expect from the cartoon of hybridizing 
the single particle excitations with excited pair states. 

Importantly, the zero temperature polarized normal 
state is a Fermi liquid. For example, the excitations near 
the Fermi surface in Fig. 2] have lifetimes which scale as 
r ~ (E — /i)~ 2 . The occupation numbers, illustrated in 
the inset of Fig . show a discontinuity. As noted by 
Haussmann [4fJ] in the early 90's, for large k, n k ~ fc~ 4 . 



This is a consequence of the discontinuity in the slope 
of the wavefunction when two atoms of unequal spin ap- 
proach one-another. More recently, Tan quantified this 
relationship, relating the coefficient of the power law to 
the "contact" 4l| . The value of this coefficient has been 
well explored both theoretically and experimentally in 
the ballanced case 42]. In the imbalanced gas, one ex- 
pects that for large k the occupation numbers should 
agree rikf = n^j.. This large k equality follows from the 
fact that the short range discontinuity involves the rel- 
ative wavefunction of an atom of each species. In the 
absence of self-consistency, Eq ([1]) fail to reproduce this 
result. However, in the theory used here (Fig. [5]), where 
we shift the chemical potentials, the tails coincide [35]. 

What are the consequences of this structure for ex- 
periment? Since the pseudogap has moved to neg- 
ative energies in the majority-species spectral den- 
sity, it can in principle be observed by final-state free 
RF-spectroscopy experiments. Unfortunately, an ex- 
plicit calculation of the momentum unresolved spectrum, 
I*(yrt) = I d 3 k/(2Tr 3 ) A ka (k 2 /2m- fi a ~ v rf )f(k 2 /2m- 
fJ-a — Vrf) does not show any particularly distinct feature 
which can be unambiguously attributed to the pseudogap 
(Fig. [S]). Momentum resolved spectra, which give infor- 
mation analogous to photo-emission [18|, directly mea- 
sure A ka {uj). As shown in Fig. [31 the finite energy pseu- 
dogap is apparent in this function as a strong energy 
dependence in the quasiparticle lifetime. 

As seen in Fig. [51 at large S the minority and ma- 
jority rf spectra coincide. This is a consequence of the 
fact that the tails of the spectrum correspond to short 
range physics, and is a measure of Tan's contact (as 
derived by several authors [HI, 0, |44| , it has the form 
I(y) = CV- 3/2 /4tt 2 \/2to). As with the tails of the oc- 
cupation numbers (Fig. [SJ, the tails of the majority and 
minority spectra fail to coincide if one does not adjust 
the chemical potentials 35]. Note that as one moves to 
the BEC side of the resonance, one experimentally sees 
that the majority species RF spectrum becomes bimodal 
45] , but near unitarity, where our calculations were per- 
formed, the experimental spectrum is unimodal. 
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